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Abstract 

The inclusive cross section for the photoproduction of neutral pions has been measured 
as a function of the transverse momentum, rapidity, and Feynman x of the ir° mesons at 
an average photon-proton centre-of-mass energy of 208 GeV and for photon virtualities 
below Q 2 = 0.01 GeV 2 . The tt° measurement extends the range covered by previous 
charged particle measurements at HERA by two units of rapidity in the photon direction 
down to a value of —5.5 in the 7p centre-of-mass frame. The 7r° transverse momentum 
distribution is well described over the whole measured range by a power law ansatz, while 
an exponential fit falls below the data at transverse momentum values above 1.5 GeV/c. 
Good agreement with the predictions of the Monte Carlo models PYTHIA and PHOJET is 
found. In the context of the PYTHIA model the data are inconsistent with large intrinsic 
transverse momentum values in the photon. 
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1 Introduction 



At the electron-proton collider HERA, photoproduction (7P) processes are induced by quasi- 
real photons that are emitted by the incoming electron or positron[|. The majority of the photon- 
proton interactions are "soft" scattering processes in which particles with limited transverse 
momenta are produced. A small fraction of yp processes, however, involve "hard" scatterings. 
In these, the photon interacts either as a single object with a parton from the proton ("direct" 
photon interactions) or via a parton-parton scattering after fluctuating into a partonic system. 
These latter "resolved" photon processes give rise to a photon remnant following the direction 
of the incoming photon. 

In the study of photoproduction scattering dynamics inclusive spectra of identified particles 
constitute an important source of information. Investigations of inclusive spectra using charged 
particles at HERA have confirmed the general features of hard photon-proton scattering out- 
lined above. In particular it was shown that in these photoproduction processes there is an 
excess of charged particles in the region of large transverse momenta, as compared to hadron- 
hadron scattering processes, and that this excess is in agreement with perturbative Quantum 
Chromodynamics (QCD) calculations, which take into account both the direct processes and 
the hard components of the resolved photon processes. These studies were carried out using the 
central tracking detectors of the experiments HI [jl], 0] and ZEUS [§]. The photon hemisphere 
was explored for charged particles down to —1.5 units of pseudorapidity 77 in the laboratory 
systemF] and with transverse momenta up to 12 GeV/c. 

The 7T° measurement described here is an extension of the previous charged particle stud- 
ies into an as yet unexplored phase space region in which the photon remnant is expected to 
dominate the particle flow. In this region charged particles could not be reliably reconstructed, 
because of limitations in acceptance of the HI tracking devices. Thus, the present n° measure- 
ment covers two additional units of rapidity y, down to y = —3.5. We note that, since at HERA 
the 7p centre-of-mass (CM) frame itself moves with an average rapidity of roughly two units 
relative to the laboratory system, this limit corresponds to y 1P ~ —5.5, as seen in the yp CM 
frame. 

The production of neutral pions has previously been extensively studied in several fixed tar- 
get experiments [Q, using both charged particle beams and photon beams, as well as in pp 
and pp collider experiments [|5p . The CM energies in these photoproduction experiments range 
up to 18 GeV. In this paper we study, for the first time at HERA, neutral pion photoproduction 
at 7p CM energies around 208 GeV. We measure the inclusive double differential n° photopro- 
duction cross section as a function of the n° transverse momentum and rapidity in the laboratory 
frame. The n° mesons are identified through their decay into two photons which are detected in 
the electromagnetic calorimeters of the HI detector. 

The paper is organized as follows. In section 2, we give a short description of the HI 
detector components that are most relevant to the 7r° cross section measurement. In section 3, we 

1 In the running periods used for this analysis, HERA was operated with a positron beam. 

2 The pseudorapidity and rapidity are defined as r\ = — lntan(0/2) and y = | In gj^ 2 , respectively, where E 
is the particle energy, p z the momentum component with respect to the z-axis and the polar angle; the z-axis is 
given by the proton direction. 
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describe the Monte Carlo models used to correct the data. Section 4 outlines the data selection 
and the analysis method, while the systematic errors are discussed in section 5. The results are 
presented and discussed in section 6. Finally, a summary is given in section 7. 



2 Detector Description 

A comprehensive description of the HI detector is given in Ref. [f7j]. The parts of the detector 
vital for this analysis are the backward calorimeter (SpaCal), the Liquid Argon calorimeter 
(LAr), the central trackers and the electron and photon detectors of the luminosity system. 
Their main properties are described here. 

The SpaCal [jSp is a lead-scintillating fibre calorimeter separated into an electromagnetic 
(EM) and a hadronic section of equal size. The EM section consists of 1 192 cells with a cross 
section of 40.5 x 40.5 mm 2 , each read out by a photomultiplier tube. The hadronic section 
consists of 128 cells with a cross section of 120 x 120 mm 2 . The EM section has a depth of 27 
radiation lengths. The polar angle range covered by the SpaCal is 153° < 9 < 178°. 

The energy resolution of the EM SpaCal section for electromagnetically interacting particles 

is 

^« *EL= e o.oio, (i) 

where E is the energy deposited in the SpaCal. The absolute energy scale, in the lower energy 
range from 0.2 GeV to 10 GeV, is known with an uncertainty of 4% 0. The resolution in 9 is 
better than 2.5 mrad for energies above 1 GeV. 

The LAr calorimeter [ ]T0[ ] is highly segmented and consists of an electromagnetic section 
with lead absorbers, corresponding to a depth varying between 20 and 30 radiation lengths, and 
a hadronic section with steel absorbers. The resolution of the LAr EM section is 

a(E) 0.12 

-^«-==©0.01, (2) 
E ^/E GeV 



as measured in test beams Q1TQ. Furthermore the absolute energy scale of the LAr calorimeter, 
in the low energy range from 0.8 GeV to 10 GeV, is known to a precision of 4%. The LAr 
calorimeter covers the polar angle range 4° < 9 < 154°. 

The HI central tracking system [0] is mounted concentrically around the beam line and 
covers polar angles in the range 20° < 9 < 160°. Measurements of the momenta of charged 
particles are provided by two coaxial cylindrical drift chambers (central jet chambers, CJC), 
mounted inside a homogeneous magnetic field of 1.15 Tesla. Further drift chambers are placed 
radially inside and outside the inner CJC chamber, providing accurate measurements of the 
z coordinates of charged tracks. Finally, multiwire proportional chambers (MWPC), which 
allow triggering on those tracks, are also located within and between the CJC chambers. In the 
present analysis tracking information is used to identify isolated electromagnetic clusters in the 
LAr calorimeter, and to find the interaction vertices of the selected events. 
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To detect photoproduction events, the electron and photon detectors of the HI luminosity 
system [[7|] are used. They consist of TICl/TlBr crystal calorimeters with an energy resolution of 
a(E)/E = 0.22/a/E/GcV and are located at z = -33 m and z = -103 m, respectively. The 
detectors serve to identify scattered positrons from photoproduction processes. The signature of 
photoproduction is a signal in the electron detector with no accompanying signal in the photon 
detector. 



3 Event Simulation 

In this analysis, photoproduction interactions are modelled by two Monte Carlo generators, 
PHOJET (version 1.04) [0 and PYTHIA (version 5.722) [jTJ]. These have both been shown to 
describe photoproduction data in previous analyses [JJ, Q [T^fl . Both generators have in common 
that they use leading order (LO) QCD matrix elements for the hard scattering subprocess. Initial 
and final state parton radiation and the LUND fragmentation model [ |T5| ] for hadronization are 



included, as implemented in the JETSET program [ jTq ] . The programs differ, however, in the 
treatment of multiple interactions and the transition from hard to soft processes at low transverse 
parton momentum. 

The PHOJET event generator simulates all components that contribute to the total photo- 
production cross section. It is based on the two-component dual parton model (for a review, see 
|Jl7p). PHOJET incorporates simulations of soft processes on the basis of Regge theory. Soft 
and hard processes are connected through a unitarization scheme which serves to keep cross 
sections finite. 

The PYTHIA event generator uses LO QCD calculations for the primary parton-parton scat- 
tering process and for multiple parton interactions. The latter are considered to arise from the 
scattering of partons from the photon and proton remnants. PYTHIA also models both hard 
and soft hadronic interactions - the latter by the exchange of low energy gluons - applying a 



unitarization scheme [18]. 



Additional transverse momentum is generated in PYTHIA by virtue of the assumed intrin- 
sic (primordial) transverse momentum (k±) of the partons in the interacting hadrons. The k± 
distribution of the quark and antiquark from the photon, fluctuating into a hadronic state, is 
parametrized as a power law function 1 / (kj_ Q + k\) [ [T8] ] with the value k±_Q = 400 MeV/c. For 
this function PYTHIA applies for the primordial k± a cutoff, the value of which is taken from 
the transverse momentum of the individual hard scattering subprocess. PYTHIA offers alterna- 
tive parametrizations for the k± distribution, by either a Gaussian or an exponential function. 
The predictions of the model for different k± distributions are compared with the experimental 
results in section 6. 

The PHOJET 1 .04 event generator has no explicit intrinsic transverse momentum for partons 
entering a hard scattering. 

For both Monte Carlo models the factorization and renormalization scales were set to the 
transverse momentum of the final state partons. GRV-LO [ [TP] ] parton distribution functions for 
the proton and photon were used for the generation of events. 
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4 Data Selection and Analysis 



Photoproduction events were selected, which have the scattered positron registered in the elec- 
tron detector of the luminosity system and consequently have a four-momentum transfer Q 2 < 
1(T 2 GeV 2 , where Q 2 = AE e E' e cos 2 (6/ '2). Here E e and E' e are the incoming and scattered 
positron energies, and 9 the polar angle of the scattered positron, 9 ~ 180°. More specifically, 
the events were required to have a positron candidate registered in the fiducial volume of the 
electron detector, with energy E' e > 4 GeV, and to have less than 2 GeV deposited in the pho- 
ton detector. The latter condition suppresses background from the proton beam appearing in 
coincidence with the high rate of events from the Bethe-Heitler reaction ep — > cyp. In order to 
ensure full efficiency and acceptance of the photoproduction tagging condition, the range of the 
variable y B = 1 — (E'J E e ) is restricted to 0.35 < y B < 0.65 in the present analysis. 

Two data samples were used in the analysis: 

• The "minimum bias" data sample, collected in 1997 and corresponding to an integrated 
luminosity of ~ 300 nb _1 . In addition to the photoproduction tagging condition, the 
trigger demanded only loose conditions on the presence of charged tracks in the central 
tracking detectors. 

• The "SpaCal" sample, collected in 1996 and corresponding to an integrated luminosity of 
~ 4.3 pb _1 . The trigger used in this sample contained, in addition to the photoproduction 
tagging condition, requirements on the energy registered in the EM SpaCal. Thus at least 
one SpaCal energy cluster (i.e. an isolated energy deposition, separated from neighbour- 
ing energy depositions by the clustering algorithm) had to exceed 2 GeV, and at the same 
time have a radial distance to the beam axis of more than 16 cm. In the subsequent anal- 
ysis a trigger efficiency of at least 40 % was ensured by demanding that in each event the 
most energetic SpaCal cluster had an energy exceeding 2.2 GeV and a radial distance to 
the beam axis exceeding 16 cm. Full trigger efficiency was reached for cluster energies 
and radial distances larger than 2.6 GeV and 22 cm, respectively. 

The "SpaCal" sample allows the study of the phase space region in which the number of 
events in the "minimum bias" sample is too small, i.e. the region where the n° has both large 
negative rapidity and large transverse momentum. 

In the further event selection the ^-coordinate z v of the event interaction vertex, recon- 
structed using the central tracking detectors, had to satisfy the condition \z v \ < 35 cm. The cut 
on z v suppresses beam-gas background events. 

The resulting "minimum bias" and "SpaCal" data samples consist of ~ 115K events and 
~ 500K events, respectively. 

In the subsequent it reconstruction, additional conditions were applied in order to accept 
energy clusters in the LAr and SpaCal EM calorimeters as photon candidates. 

• For a LAr cluster to be accepted as a single photon candidate, the transverse radius of the 
cluster had to be smaller than 8 cm and the longitudinal cluster extension below 10 cm. 
The cluster energy had to exceed 0.3 GeV. Furthermore, the distance between the cluster 
position and the impact point of the closest reconstructed charged track had to be larger 
than 5 cm. 
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• An EM SpaCal cluster was associated with a photon if it was reconstructed as an isolated 
energy deposition covering at least two SpaCal cells. The cluster energy had to exceed 
0.3 GeV. A radial distance between 8 cm and 75 cm from the cluster to the beam axis 
was required to ensure full energy containment. Information from the backward tracking 
detectors was not used in selecting photon candidates in the SpaCal, since many photons 
shower in the material in front of the backward tracking detectors and the SpaCal. 

The procedure for reconstructing and counting tt° mesons was as follows: Each photon pair, 
from either the LAr or the SpaCal and satisfying the above criteria, was labelled according to 
transverse momentum p±_ and rapidity y. Photon four-momenta were derived from the cluster 
energy, the cluster coordinates and the coordinates of the interaction vertex. For every photon 
pair, the invariant mass m 77 was calculated according to 

m 77 c 4 = 2E 1 E 2 (1 - cos ^i 2 ) . (3) 

Here E\ and E 2 are the energies of the measured calorimeter clusters, attributed to the decay 
photons, and is their opening angle. Fig. |T] shows two-photon invariant mass distributions 
in several regions of transverse momentum p± and rapidity y of photon pairs measured in the 
SpaCal (Fig. [l]a and [l]b), and in the LAr calorimeter (Fig. [T]c and [TJd). These variables were 
calculated from the four-momentum sum of the two photons. 

In each bin of p± and y the two-photon mass distribution can be satisfactorily described 
by the sum of a Gaussian function for the ir° signal and a background curve. The background 
was determined by fitting a polynomial of fourth order multiplied by a function of the type 
(m 77 — m 77 r ) a in order to improve the threshold description where the background rises steeply. 
m 77 r is given by the centre of the lowest non-empty interval in the mass distribution, and a is a 
fitted parameter with < a < 1. 

The double differential cross section for inclusive n° photoproduction is given by the ex- 
pression 

d 2 (T lp = N piod (Ap ± , Ay) 
dp\dy 2p ± ■ Ap L ■ Ay • $ ■ £ ' 

-Wprod is me number of 7r°s produced in the transverse momentum (Ap±) and rapidity (Ay) in- 
terval considered, determined with all efficiency and acceptance corrections. The efficiency was 
determined by simulation of photoproduction events using the PHOJET and PYTHIA genera- 
tors (see section 3), with full simulation of the detector response. £ is the integrated luminosity, 
and $ = 0.00968 is the photon flux factor, calculated according to the Weizsacker - Williams 
formula [|2(]] for the y B range covered in this analysis. This range corresponds to the ^p CM 
energy range 177 < y/s^ < 242 GeV, with an average of 208 GeV. We note that in the various 
bins of p± and y the average ^fs^, may deviate from this value by up to ±3%. 

5 Systematic Errors 



The errors on the results are dominated by systematic effects with the following main sources: 
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1) The uncertainty of the fitting procedure arising from the determination of the combinatorial 
background in the fits of the two-photon mass distributions. In all kinematical intervals relevant 
for this analysis this was estimated by fitting these distributions, varying the fit parameters 
and fit ranges several times, following a random choice. Also the order of the polynomial 
used to describe the background was varied. The errors derived in this way vary substantially 
due to the statistical uncertainty of the background estimation. In order to get a conservative 
but reliable estimate for the uncertainty from the fitting procedure, a global error of 8% was 
therefore assigned, based on the distribution of the errors. 

2) The 7t° reconstruction efficiencies, derived using the PYTHIA and PHOJET models, are in 
general in good agreement, except in the bins of high p± where the two models differ by up to 
30%. Such differences can be understood as being due to internal differences of the models in 
the population of certain kinematical regions. In the cross section calculations we have taken 
the mean reconstruction efficiency and assigned the difference as a bin to bin systematic error, 
varying between 1% and 15%. 

3) The uncertainty of the energy scales for the LAr and SpaCal calorimeters: Shifting the scale 
upward (downward) by 4% increases (reduces) the cross sections by about 20% for the data 
selected with the SpaCal trigger, and by about 10% for the minimum bias triggered data. 

4) The uncertainty on the electron detector acceptance is 4% and 6%, for the 1996 and 1997 
data samples, respectively. 

5) The uncertainty on the efficiency of the SpaCal trigger is found to be 1%, from cross checks 
with an independent trigger. 

6) The contribution of 7r°s from beam-gas interactions has been estimated to be less than 1% 
using the so-called "pilot bunches", i.e. positron beam bunches without partners with which 
they can collide in the proton beam. 

7) The integrated luminosity is known to a precision of better than 2%. 



6 Results 

The cross section for inclusive neutral pion production in yp interactions is displayed in Fig. [2| 
as a function of the n° rapidity for four intervals of the transverse momentum and in Fig. |3] as 
a function of p± for five intervals of rapidity. The data are summarized in Table 1 . Rapidity 
is measured in the laboratory frame. Bin-centre corrections have been applied and errors are 
quoted as the sum of statistical and systematic errors, added in quadrature. In regions of low 
statistics (low y, high p±) or large combinatorial background (central region, low p±) unstable 
fits prevent the determination of cross section values. 

Values derived from the HI charged particle cross section measurements [[]]], which cover 
the pseudorapidity range |^| < 1.5, are also shown (full triangles) in Figs. ^ and ^ In order to 
compare these measurements with the n° result, the charged particle cross sections were first 
corrected for the non-pion part by subtracting a fraction of (17.5 ± 0.2)% as derived from the 
PYTHIA and PHOJET predictions (the error reflects the difference in the predictions). The 
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resulting n + + 7r~ cross sections were then further divided by the isospin factor 2. The ep 
cross sections given in [[]]] have been converted into ^p cross sections using the appropriate flux 
factor. While the charged pion cross section entries in Fig. ^correspond directly to the original 
data points in [jl]], a smooth interpolation was used to arrive at the charged pion cross sections 
in Fig. ^, in order to account for the different transverse momentum ranges. The agreement 
between the neutral and charged pion cross sections is good in the kinematical range covered 
by both analyses. Figs. ^ and ^] demonstrate well the additional phase space region accessed 
with the present neutral pion analysis. 

The model predictions according to the PHOJET and PYTHIA simulations are also shown 
in Figs. |2| and|3], with satisfactory agreement especially for PYTHIA. PHOJET predicts slightly 
too large a n° rate at high p±_ and low y. 

The 7T° transverse momentum spectrum, integrated over the laboratory rapidity range cov- 
ered by the SpaCal, i.e. —3.5 < y < —1.5, is shown in Fig. 0. The data are summarized in 
Table 2. In deriving these cross sections also those regions were included for which no results 
are given in Table 1. 

Transverse momentum spectra in high energy hadron - hadron collisions are successfully 
described by a power law ansatz of the form 



This QCD inspired ansatz pl| , [221 ] was designed to describe transverse momentum spectra of 
centrally produced particles. It was shown in [fl], 0, ^ that it fits well the HI and ZEUS spectra 
measured for charged particles in the central pseudorapidity range. This is in agreement with 
observations made in many hadron-hadron experiments, see e.g. Kourkoumelis et al. ^\ for 
centrally produced 7r°s at various centre-of-mass energies at the CERN ISR, and Banner et 
al. [Bp for measurements performed at ^Js = 540 GeV at the CERN SPS. A review of large pj_ 



particle production at the ISR is given by Geist et al. [23]. 



y 


p ± [GeV/c] 


[0.2,0.6] 
(0.36) 


[0.6,0.8] 
(0.69) 


[0.8,1.0] 
(0.89) 


[1.0,2.0] 
(1.31) 


[-3.5,-2.8] (-3.15) 


47 ± 10 


5.0 ± 1.2 






[-2.8,-2.4] (-2.6) 


94 ± 16 


12 ±2 


3.6 ±0.8 


0.23 ±0.08 


[-2.4,-2.0] (-2.2) 


158 ± 27 


27 ±5 


8.0 ±2.0 


0.68 ±0.17 


1-2.0,-1.5] (-1.75) 


176 ± 27 


28 ±5 


12.9 ±3.4 


1.34 ±0.32 


[-0.5, +1.0] (0.25) 






16.9 ±3.5 


1.97 ±0.35 



Table 1: Inclusive n° photoproduction cross sections d 2 a 1P / ' dp\dy (yub/(GeV/c) 2 ) for bins in 
the 7T° transverse momentum p± and rapidity y. The bin centres are given in addition to the 
interval limits. 
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y 


p± [GeV/c] 


[0.2,0.6] 
(0.34) 


[0.6,0.9] 
(0.73) 


[0.9,1.2] 
(1.03) 


[1.2,1.5] 
(1.33) 


[1.5,1.8] 
(1.63) 


[1.8,2.1] 
(1-93) 


[-3.5,-1.5] 


114 ±20 


12.5 ±2.4 


2.0 ±0.4 


0.47 ±0.15 


0.14 ±0.05 


0.04 ±0.01 



Table 2: Inclusive ir° photoproduction cross sections d 2 a yp j dp\dy (/ib/(GeV/c) 2 ) for bins in 
the 7T° transverse momentum p± and for the rapidity range —3.5 < y < —1.5. The bin centres 
are given in addition to the interval limits. 



The measurements presented in Fig. |] have been performed in the backward region at large 
negative rapidities where phase space effects begin to be visible, causing a damping of trans- 
verse momentum spectra. This leads to a reduction of the cross section and to a steepening of 
the p± distribution. The power law ansatz does, nevertheless, satisfactorily describe the shape of 
the spectrum if the parameter p± is fixed to the value found in [|1]], namely p± = 0.63 GeV/c. 
A fit with the power law ansatz yields n = 8.0 ± 0.2 (x 2 /ndf = 0.85). The value is slightly 
larger than the values found in the fits to the charged particle spectra in the central rapidity re- 
gion, n = 7.1 ± 0.2 [[Ufl n = 7.03 ± 0.07 0] (for p x > 2 GeV/c) and n = 7.25 ± 0.03 [|] (for 
p± > 1.2 GeV/c, withp ±0 = 0.54 GeV/c). 

The measured p±_ distribution in Fig. |3] can also be described by an exponential distribution 
of the form 



d 2 a 
dp\dy 



a ■ exp(— b ■ \ / p 2 ± c 2 + m 2 c A ). (6) 



The fit to the data in the range 0.3 < p± < 1.6 GeV/c is shown in Fig. 4a and yields the 
value b = (5.7 ± 0.3) GeV 1 (% 2 /ndf = 1.0). Such slope values are typical of soft hadronic 
interactions and expected in thermodynamical models (see e.g. [p2|]). Above p± « 1.5 GeV/c, 
however, the exponential form tends to fall below the data, in contrast to the power law ansatz. 
This observation can be interpreted as an indication of the onset of hard parton-parton scattering 
processes which lead to a hardening of the transverse momentum distribution. 

This interpretation is supported by a comparison of the data with the MC models. Fig. ^b 
shows again the cross section as a function of p±, this time together with the model predictions 
of PYTHIA and PHOJET. The agreement between the data and the PYTHIA prediction is in 
general slightly better than for PHOJET. Also shown is the PYTHIA prediction for the amount 
of "direct" photon interactions. The contribution of this direct component accounts for ~ 5% 
of the cross section at the lowest p± values, increasing up to ~ 20% at the highest p± values 
accessed. Thus, in the rapidity and transverse momentum range covered by these measurements, 
the cross section is dominated by "resolved" photon processes. 

The generator PYTHIA has also been used to search for evidence for a non-zero intrinsic 
transverse momentum, k±, of partons inside the photon. As mentioned in section 3, PYTHIA 
offers the choice of several distributions for the parameter k±: a Gaussian, an exponential, and 
a power law ansatz. All three parametrizations have been used and the parameter k± , which 



3 In Ref. |JtJ] the error on the power n is erroneously quoted as 2.0. 
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defines their widths, varied between and 2 GeV/c. The largest influence on the predicted shape 
of the y distribution appears in the highest pj_ interval. This interval, together with the curves 
predicted by PYTHIA for the Gaussian, exponential, and power law parametrizations, is shown 
in Fig. |5|. We note that for these plots no additional cutoff for the primordial k± has been applied 
(see section 3). For all three functions large effects are seen when varying the parameter k±o, 
and low values of k± are preferred. Even a vanishing k± is consistent with the data. 

Finally the behaviour of the distribution of the Feynman variable xp has been investigated, 
with xp defined as the fraction pl/Pl 4 ""* where pl is the pion's momentum component parallel 
to the beam axis in the jp CM frame, and p™ ax its maximum value in a given event, namely 

n max ^ rr, /Q 

Pl ~ v 7P' 

In Fig. the differential cross section da lp /dxp is displayed, obtained in the rapidity range 
—3.5 < y < —1.5. There is good agreement with the model predictions of PYTHIA and 
PHOJET. The cross section values are given in Table ^ 



y 


Xp 


[0.1,0.15] 
(0.125) 


[0.15,0.2] 
(0.175) 


[0.2,0.25] 
(0.225) 


[0.25,0.3] 
(0.275) 


[0.3,0.4] 
(0.348) 


[0.4,0.6] 
(0.493) 


[-3.5,-1.5] 


477 ± 67 


279 ± 43 


209 ± 35 


131 ± 19 


65 ± 12 


19 ±5 



Table 3: Inclusive n° photoproduction cross sections da 1P /dxp (pb) in intervals of Xp, in the 
rapidity range —3.5 < y < —1.5. The bin centres are given in addition to the interval limits. 



Data from the Omega Photon Collaboration (WA69, Apsimon et al., [|5|]) are also shown 
in Fig. |[ They were obtained in fixed target jp collisions with an average photon energy of 
80 GeV (corresponding to ^fs^, = 12.3 GeV). From this CM energy and the 7r° mass one 
derives a CM rapidity range of up to 4.5 units available for neutral pions. These cross sections 
of inclusive ir° photoproduction, available in Ref. [ |2"4"| ] as Ed 3 a/dp 3 in bins of xp and p±, have 
here been converted to cross sections differential in xp, by integrating over p± in each brnP] of 
Xp. The two experiments are in good agreement, although small differences are seen in shape 
and normalization. These differences can be qualitatively understood from the difference in 
rapidity range; the cuts —3.5 < y < —1.5 (corresponding to — 5.5< y < — 3.5 in the ^p CM) 
reduce the HI n° cross section at both large and small xp, a reduction which is only partly 
compensated for by the expected increase due to the rise of the total cross section between the 
two CM energies (factor ~ 1.3 p5|]). 

7 Summary and Conclusions 

Inclusive n° photoproduction in the photon hemisphere has been analysed at HERA in the 7p 
CM energy range 177 < y^s^ < 242 GeV, with average ~ 208 GeV. Differential cross sections 
with respect to transverse momentum p±, rapidity y and the Feynman variable xp are presented. 

4 In the lowest bin of p± and xp, where a measurement with 80 GeV beam energy is missing in [^[j, the 
corresponding measurement from the data set with average beam energy 140 GeV has been used. 
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The neutral pion data extend the measurements towards the previously unexplored domain of 
large rapidity in the photon fragmentation region, dominated by resolved -yp interactions. This 
work extends previous HI measurements performed with charged particles. In the phase space 
region covered by both analyses there is good agreement between the neutral and charged pion 
differential cross sections. 

The differential cross section as a function of transverse momentum p ± shows an exponen- 
tial fall at lower pj_ values as seen in soft hadron-hadron collisions, but exhibits at values of 
p± larger than 1 .5 GeV/c an enhancement which is expected for hard parton-parton scattering 
processes. The distribution in the entire p± range covered here is well described by a power law 
ansatz. The event generators PHOJET and PYTHIA are able to describe the measured cross 
sections in this kinematical domain, with a slight preference for PYTHIA. In the context of the 
PYTHIA model, the data are inconsistent with large values of an intrinsic transverse momentum 
in the photon. 
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Figure 1 : Two-photon invariant mass distributions reconstructed in the SpaCal calorimeter (a, 
b) and in the Liquid Argon calorimeter (c, d). The full curves show fits of a sum of a Gaussian 
and a background distribution (dashed curves) as described in the text. 
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Figure 2: Inclusive tt° photoproduction cross section as a function of laboratory rapidity y in 
intervals of transverse momentum p± (full circles). The triangles are the corresponding cross 
section values for charged pions, derived from [|T|] by subtracting a fraction of 17.5% to account 
for the "non-pion" contribution, and then dividing by the isospin factor 2 (see text). The curves 
are the predictions of the PHOJET (dashed) and PYTHIA (full) event generators for neutral 
pions. 
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Figure 3: Inclusive n° photoproduction cross section as a function of transverse momentum 
p±_ in intervals of laboratory rapidity y (full circles). The triangles are the corresponding cross 
section values for charged pions in the pseudorapidity range \r}\ < 1.5, derived from [JTJ] by 
subtracting a fraction of 17.5% to account for the "non-pion" contribution, and then dividing 
by the isospin factor 2 (see text). The curves are the predictions of the PHOJET (dashed) and 
PYTHIA (full) event generators for neutral pions. 
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Figure 4: (a) Inclusive n° photoproduction cross section as a function of transverse momentum 
p± (circles), for the rapidity region —3.5 < y < —1.5. The dash - dotted curve is an exponential 
fit to the data in the range 0.3 < p± < 1.6 GeV/c; the full curve is the result of a power law fit to 
all data points (see text), (b) Same cross section as in (a) compared to predictions of PHOJET 
and PYTHIA. In addition, the contribution of direct photon interactions in PYTHIA is shown. 
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Figure 5: The n° cross section as a function of rapidity y for the transverse momentum range 
1.0 < p± < 2.0 GeV/c. The curves are the predictions of PYTHIA, labelled with the corre- 
sponding k± values, as described in the text. No additional cutoff on k ± is applied (see text). 
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Figure 6: Inclusive n° photoproduction cross section as a function of Feynman x. The curves 
are the predictions of the PHOJET (dashed) and PYTHIA (full) models. The HI data, at ^/s^, ~ 
208 GeV, as well as the MC model predictions, are obtained in the rapidity range —3.5 < y < 
— 1.5. The data of the Omega Photon Collaboration (WA69), at ^fs^, m 12.3 GeV, were derived 
as described in the text. 
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